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We investigate the role of disorder on the superconducting (SC) fluctuation in short coher- 
ence length d-wave superconductors. The particular intetest is focused on the disorder-induced 
microscopic inhomogeneity of SC fluctuation and its effect on the pseudogap phenomena. We 
formulate the self-consistent 1-loop order theory for the SC fluctuation in inhomogeneous sys- 
tems and analyze the disordered t-t'-V model. The SC correlation function, electronic DOS and 
the critical temperature are estimated. The SC fluctuation is localized like a nanoscale gran- 
ular structure when the coherence length is short, namely the transition temperature is high. 
This is contrasted to the long coherence length superconductors where the order parameter is 
almost uniform in the microscopic scale. In the former case, the SC fluctuation is enhanced 
by the disorder in contrast to the Abrikosov-Gorkov theory. These results are consistent with 
the STM, NMR and transport measurements in high-Tc cuprates and illuminate the essential 
role of the microscopic inhomogeneity. We calculate the spacial dependence of DOS around the 
single impurity and discuss the consistency with the NMR measurements. 

KEYWORDS: Microscopic inhomogeneity; pseudogap; superconducting fluctuation; high-Tc cuprates; ran- 
domness 



1. Introduction 

The pseudogap phenomena in high-Tc cuprates have 
been one of the central issues in strongly correlated elec- 
tron systems. ""^ The theoretical elucidation of anomalous 
properties is highly desired for a comprehensive under- 
standing of high-Tc superconductivity. 

Among many theoretical proposals, the "pairing sce- 
nario" is one of the promising ideas. Then, the origin 
of pseudogap is attributed to the superconducting cor- 
relation above Tc.^"^'* In contrast to the conventional 
superconductors, the coherence length is extraordinary 
short in under-doped cuprates, which is typically a few 
lattice spacing. The superconducting (SC) fluctuation is 
enhanced in such a short coherence length superconduc- 
tor with quasi-2D lattice structure. 

In order to examine the validity of pairing scenario it 
is important to investigate the roles of disorder and/or 
magnetic field which suppress the d-wave superconduc- 
tivity. As for the magnetic field, it has been shown that 
the pseudogap induced by the SC fluctuation is slightly 
affected by the orbital de-pairing effect around the on- 
set temperature T = T*.^^ This is mainly because the 
correlation length of superconductivity is still short at 
T = T*, where the pseudogap is induced by the short 
range SC correlation. On the other hand, the critical 
fluctuation around T = Tc is remarkably affected by the 
magnetic field because the correlation length diverges at 
the critical point. ^^'^^ These results are consistent with 
NMR measurements.^^' 

In this paper, we investigate the role of disorder on the 
short coherence length superconductor with cZ-wave sym- 
metry. This issue is particularly interesting because the 
microscopic inhomogeneity observed in recent STM mea- 
surements^^"^^ indicates a novel effect of disorder in the 



under-doped region. Furthermore, it is considered that 
the disorder plays an essential role in the superconductor- 
insulator (SI) transition in high-Tc cuprates as implied 
by the disordered magnetism in the LSCO and YBCO 
compounds. ■^^"^^ 

The pseudogap induced by the d-wave SC fluctuation 
has been investigated in disordered system on the basis of 
the self-consistent t- matrix approximation (SCTMA)^^ 
for the disorder. Then, it was shown that the pseudo- 
gap and SC fluctuation are significantly suppressedhy the 
disorder. This conclusion seems to be reasonable since 
the Tc in non-s-wave superconductors is suppressed by 
the disorder. Then, the thermal fluctuation can be sup- 
pressed. However, this result is incompatible with the 
experimental results which have indicated a robustness 
of pseudogap in disorder-doped materials. '^""'^■^ 

In this paper, we carry out the calculation beyond the 
SCTMA where the disorder is exactly taken into account. 
Contrary to the SCTMA, it is shown that the SC fluc- 
tuation is enhanced by the microscopic inhomogeneity 
which is characteristic in the quasi-2D short coherence 
length superconductor with non-s-wave symmetry. Then, 
the long range coherence hardly develops owing to the 
nanoscale inhomogeneity of SC order parameter while 
the well developed short range correlation leads to the 
pseudogap in the low energy spectrum. This result illu- 
minates the breakdown of Abrikosov-Gorkov theory for 
the short coherence length superconductors. 

We analyze the t-t'-V model to investigate the mi- 
croscopic role of disorders on the SC fluctuation. We 
adopt the I-loop order approximation with respect to 
the SC fluctuation, which is called "self-consistent T- 
matrix approximation" This approach is comple- 
mentary with the phase fluctuation theory^^^"^ which is 
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a phenomcnological description for the deeply critical re- 
gion. It has been shown that the 1-loop order theory is 
consistent with many experimental results in the pseu- 
dogap state. The anomalous behaviors in single parti- 
cle, magnetic and transport properties^' are explained 
in a comprehensive way by taking into account the spin 
fluctuation in addition to the SC fiiictuation.^'*''^^''^^ Fur- 
thermore, the microscopic estimation of the SC fluctua- 
tion in Hubbard model is consistent with the remarkable 
electron-hole asymmetry in high-Tc cuprates.^^ 

This paper is organized as follows. The disordered t- 
t'-V model is introduced in §2, where we formulate the 
self-consistent T-matrix approximation^^ in disordered 
systems. In §3.1, the typical results on the spatial depen- 
dence of SC fluctuation are shown. In §3.2, we numeri- 
cally take the random average and calculate the DOS, 
SC correlation function and Tc- The role of SC fluctua- 
tion is clarified by comparing with the mean field theory. 
The effects of microscopic inhomogeneity are illuminated 
by the comparison with the SCTMA for disorder effects. 
In §4, we briefly discuss the spatial dependence of DOS 
around the single impurity. Some discussions are given 
in the last section §5. 

2. Formulation in the Disordered t-t'-V Model 

2.1 1-loop order theory in disordered systems 

We adopt the disordered t-t'-V model which is ex- 
pressed as, 

H = -t c\^Cja+t' 

<i,j>,a <Si,jS>,CT 

<i,j> i 

where is the electron number at the site i. The sym- 
bols < i,j > and <^ i,j denote the summation over 
the nearest neighbor sites and that over the next near- 
est neighbor sites, respectively. The attractive interac- 
tion V < between the nearest neighbor sites favors 
the (ix2-y2-wave superconductivity around the half fill- 
ing. The disorder is introduced by the disordered poten- 
tial Wj. Although we have investigated several kinds of 
disorder including the box disorder, Gaussian disorder 
and t'-disorder, we show the results for the point disor- 
der where Wi = or Wi = W. Since many experimental 
studies on the Zn-doped cuprates have been reported, 
we assume the strong disorder W = AOt ^ t. Then, the 
sites i where Wi = W correspond to the Zn sites and the 
others are the Cu sites. Recently, the disorder outside 
the plane has attracted much interests.^^'^^'^^ Then, the 
disordered potential is described by the extended point 
disorder. We have confirmed that the following results 
are qualitatively independent of the type of disorder. 

It should be noted that the t-t'-V model is a semi- 
phenomenological model for the short coherence length 
d-wave superconductor. The d-wave pairing in high-Tc 
cuprates is originally induced by the short range re- 
pulsive interaction through the many body effect. We 
have analyzed the repulsive Hubbard model within the 
FLEX-I- T-matrix approximation,^^ and shown that the 
attractive models like eq. (1) appropriately capture the 
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role of SC fiuct nation.^"' In this paper, the effect of the 
microscopic inhomogeneity will be elucidated from the 
general point of view. 

We choose the unit of energy as t = 1 and fix t' /t = 
0.25. The chemical potential fi is fixed to ^ = —0.8 which 
corresponds to the 10% hole doping. The typical Fermi 
surface of under-doped cuprates is well reproduced by 
this parameter set. Since we fix the chemical potential 
instead of the total electron number, the average electron 
density per Cu site n = X]<i> '^i/-^o is little affected by 
the disorder-doping. Here, denotes the summation 

over the sites where Wi = Q. Nq ]s the number of those 
sites. Although the effect of impurity doping on the aver- 
age electron density in high-Tc cuprates is not clear, we 
can study the general roles of disorder by analyzing this 
model. 

We define the non-interacting part of Hamiltonian as 
follows, 

<i,j>,<7 «i,j»,<7 i 

= ^ct^0Sa. (2) 

Here, cj. — (c|g.) is the TV-component vector and N = 
T X L is the total number of sites. The Green function 
in the non-interacting system {V = 0) is obtained by 
diagonalizing the N x N Hermitian matrix Hq as, 

Go{i,j,uin) = y^^u{m,l)- — u*{n,l), (3) 

ILOn - ei 

where ei is the l-th eigenvalue of the matrix Hq and 
u{m, I) is the m-th element of the eigenvector. The in- 
dex m (n) is defined so that the m-th (n-th) element 
of cl is the creation operator at site i (j). The fermion 
Matsubara frequency is given as uin = (2n-|- l)7rT. Since 
we have included the disorder potential into the unper- 
turbed Hamiltonian, the disorder is exactly taken into 
account in the following calculation. 

In ease oi V ^ 0, the dressed Green function is ob- 
tained by the Dyson equation as, 

G(l,J,Un) = Go{i,j,UJn) + 

YGo{i,i',u>n)^(i',j',u;n)G{j',j,u;n), (4) 

where I](z,j,a;„) is the self-energy. In this paper we es- 
timate the self-energy within the self-consistent or non- 
sclf-consistent T-matrix approximation whose results in 
the clean system have been summarized in Ref. 14. 

It generally takes a long computational time to calcu- 
late the self-energy in the inhomogeneous system because 
the Green function G{i,j,uin) depends not only on the 
relative coordinate i — j but also on the center of mass 
coordinate [i + j)/2. In order to reduce the computa- 
tional time, we apply the quasi-static approximation for 
the propagator of SC fluctuation, namely the T-matrix. 
The quantum dynamics of SC fluctuation, namely the 
imaginary time dependence of fluctuation propagator, is 
ignored in this approximation and only the thermal fluc- 
tuation is taken into account. ^^"^^ This approximation 
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Fig. 1. The diagrammatic representation of (a) T-matrix and (b) 
self-energy. 



is valid at finite temperature around the critical point. 
We have used the quasi-static approximation for the es- 
timation of higher order corrections beyond the 1-loop 
order.*® Then, we have found that this approximation is 
quantitatively justified in the pseudogap region. It should 
be noted that the quasi-static approximation is partic- 
ularly appropriate for the SC fluctuation owing to the 
characteristic dynamical property.*® 

Since the T-matrix having non-zero Matsubara fre- 
quency 0„ 7^ is ignored in the quasi-static approxi- 
mation, we simply drop the variable of frequency in the 
T-matrix. Then, the T-matrix is obtained by the follow- 
ing equations, 



T{i,j,5,5') = To{i,j,Sj')- 

To{lk,S,S")VTik,j,5",S'), 



k,5" 



(5) 



where 5 = (±1,0) or j = (0, ±1). The diagrammatic 
representation of T-matrix is shown in Fig. 1(a). 

The self-energy is represented in Fig. 1(b) and the ex- 
pression is given as. 



S(U-,a;„) = -y2^[T,(i,j,5,5') +3Tt(i,j,5,5') 

S,5' 

xG{j + 5' ,i + 5,-ojn), 



(7) 



where 

T,,,(i,j,l5') = \t{1115') ± ^T(i;j + 5',5,-5'). (8) 

Here, Ts{i,j,S,S') (Tt{i,j,S,5')) is the fluctuation prop- 
agator for the spin singlet (triplet) pairing. The former 
includes the rfx2_y2-wave and extended s-wave pairings 
while the latter includes the p^-wavc and py-wave pair- 
ings. The coeflacient 3 of Ti{i,j, 6, 6') in eq. (7) is due to 
the spin 1 of Cooper pairs in the spin triplet pairing. 

Although the d-wavc pairing state is most stable in 
the t-t'-V model, the pairing correlation develops also 
in the p-wave channel. The latter is expected to be an 
artifact of the t-t'-V model because the p-wave pairing 
instability is negligible in the microscopic model, such as 
the Hubbard, d-p and t-J models.^^'*'''*^ Therefore, we 



simply drop the contribution from the spin triplet SC 
fluctuation in eq. (7) and obtain the self-energy as, 

E(i, 1 cOn) = -V^ J2 ^^(^' ^')G(i + + 6, -w„).(9) 

6,5' 

It should be again stressed that the t-t'-V model is a 
semi-phenomenological model which is relevant for the 
d-wave SC fluctuation. Note that the extended s-wave 
pairing fluctuation is also included in Ts{i, j, S, 6') and we 
cannot separate these pairings in the disordered system 
owing to the violation of 4-fold rotational symmetry. We 
have confirmed that the contribution from the extended 
s-wave pairing is not important. 

In the self-consistent T-matrix approximation, the self- 
energy, Green function and T-matrix are obtained by 
solving the self-consistent equations, namely eqs. (3-6) 
and (8-9). The non-self-consistent T-matrix approxima- 
tion is given by replacing the Green function G in eqs. (6) 
and (9) with the unperturbed one Gq. The self-consistent 
calculation is carried out for up to A'' = 31 x 31 sites, 
while the non-self-consistent one is performed for up to 
TV = 65 X 65 sites. 

2.2 SCTMA for disorder 

The superconductivity in the disordered alloy has been 
investigated for more than four decades.*^ The Born 
approximation or SCTMA have been used to discuss 
the macroscopic properties. Then, the random average 
is taken in advance and the translational symmetry is 
reestablished in the results. The role of strong disorder 
on the gap-less superconductivity has been elucidated 
with use of the SCTMA.^°"^^ However, the role of inho- 
mogeneity is not appropriately taken into account in the 
SCTMA because the interference effect of many disorders 
is neglected. In contrast to that, the disorder is exactly 
taken into account in the formulation in §2.1. Therefore, 
it is interesting to illuminate the role of microscopic in- 
homogeneity by the comparison with SCTMA. In this 
subsection, we briefiy explain the formulation of SCTMA 
in the t-t'-V model. 

Since the disorder is taken into account in an approx- 
imate way, the unperturbed Hamiltonian is defined as, 

H^ = -t Yl 4aCja + t' Yl 4cia-M$^ni. (10) 

<i,j>,(7 <<i,j>>,a i.a 

Then, the Green function, self-energy and T-matrix can 
be written in the fc-space representation. For example, 
the Green function is written as,^^' 

G{k,UJr:) = -. „ ], . ^ t y (11) 

\ujn - e{k) - Ss(fe,a;„) - Sd(w„) 

where £(fe) is the dispersion relation, 

e{k) = — 2t(cos fcx + cos ky) + At' cos fcx cos ky — fi. (12) 

The self-energy is given by the c;ontribution from the 
SC fluctuation Es(fe,a;„) and that from the disorder 
Ed(w„). In order to compare with the formulation in 
§2.1, we calculate the self-energy Ss(fc,a;„) in the same 
approximation as in §2.1. Then, T,s{k,LOn) is expressed 
in the real space as, 

S,(i, LOn) = -V^ J2 -6 + 6', -LOn), (13) 

5,6' 
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where 



TS, S, 6') = ^T{1 5, S') + ^T{i - 6', 6, -6% (14) 



and 



T{i,S,S')=ToilS,6')- 

^To{i-k,6,S")VTCk,S",5'), (15) 



k,6" 



Toil 6, S') = tY^ Gil u>m)G{i + 5-6', -a;„). (16) 

m 

The fc-space representations of T,s{k,uin) and G{k,u)n) 
are obtained from the Fourier transformation of Sg(i, a;„) 
and G{i,uin), respectively. 



// W 
// \\ 



Fig. 2. Self-energy arising from the impurity scattering in the 
SCTMA. 



The self-energy arising from the disorder is given as, 

where riimp = {N — Nq)/N is the impurity concentration. 
The diagrammatic representation of cq. (17) is shown in 
Fig. 2. Here, we have taken the limit — > oo for sim- 
plicity, which is consistent with the sufficiently strong 
potential W = 40i in §2.1. We have ignored the vertex 
correction arising from the disorder on the T-matrix be- 
cause it is negligible for the rf-wave SC fluctuation.^^' 

If wo neglect the self-energy due to the SC fluctuation 
T,s{k,uin), this calculation is reduced to the SCTMA in 
the mean field theory. The deviation from the pure d- 
wave superconductivity in the experimental results, such 
as the NMR 1/TiT and magnetic field penetration depth, 
has been resolved within the mean field theory. ^'^"^^ The 
consistency between the theory and experiment has been 
regarded as an evidence for the d-wave pairing in high-Tc 
cuprates.^^'^^ We show that this approximation breaks 
down in the under-doped region, while the SCTMA will 
be valid in the over-doped and electron-doped region. 

3. Microscopic Inhomogeneity and SC Fluctua- 
tion 

3.1 Microscopic inhomogeneity 

First of all, we show the typical spatial dependence 
of SC order parameter in the disordered short coherence 
length superconductor. In order to discuss the long range 
correlation, we show the results of non-sclf-consistent T- 
matrix approximation for = 65 x 65 sites in this sub- 
section. We have confirmed that qualitatively the same 
results are obtained in the self-consistent calculation. 

We determine the spatial dependence of SC order pa- 
rameter from the fluctuation propagator, namely the T- 
matrix. Although the amplitude of SC order parameter 



is not determined by the T-matrix, the spatial depen- 
dence is determined in the following way. First, we solve 
the eigenvalue equation to diagonalize the irreducible T- 
matrix as. 



tnv{n,i,5) = To{i,j,6,5')v{n,j,6'). 



(18) 



Here, i„ is the eigenvalue of T^ii, j ,5,5') and we define 

ti > t2 > ts > t4N- Then, the T-matrix is described 

as, 

T{i, J, 5, 6') = ^ v{n, 1, 6)^^^v{n, j, 5'), (19) 

and the critical temperature is determined by the crite- 
rion \V\ti = 1. The order parameter of d-wave supercon- 
ductivity just below Tc is obtained as. 



Ai(i) = J2i-^fva,l5), 



(20) 



where P = (P = 1) for 5 = (±1,0) (5 = (0,±1)). 
The extended s-wave component of the order parameter 
is also determined as. 



As(z)=^^;(l,z,J). 



(21) 



Although As(«) is finite in the disordered system, it is 
much smaller than the d-wave component Ad(i). The 
small s-wave component does not play any important 
role in the following results. 

We show the typical spatial dependence of Ad(i) in 
Fig. 3. We choose the parameter asV = —0.6 in Fig. 3(a) 
while we choose V = —1.5 in Figs. 3(b) and (c). The 
disorder concentration is chosen to be nimp = 0.01 in 
Figs. 3(a) and (b) and nimp = 0.05 in Fig. 3(c). Note that 
the coherence length in the clean limit is scaled by 
as ^0 vf/T^ where T° is the transition temperature in 
the clean limit and vf is the Fermi velocity. We obtain 
T° = 0.0418 and = 0.269 for V = -0.6 and V = -1.5, 
respectively. Thus, the coherence length in Fig. 3(a) is 
much longer than that in Figs. 3(b) and (c). 

We clearly see the remarkable spatial inhomogeneity 
of SC order parameter in Figs. 3(b) and (c). In partic- 
ular, the SC order parameter is highly localized by the 
5% disorder-doping in the case of V = —1.5. It should be 
noticed that these results are very unusual. In case of the 
conventional long coherence length superconductor, the 
SC order parameter is almost uniform in the microscopic 
scale as shown in Fig. 3(a). If we increase the disorder 
with fixing V = 0.6, the superconductivity is depressed 
before the SC order parameter is localized. This is the 
usual behavior of non-s-wave superconductor where the 
superconductivity is suppressed by the non- magnetic im- 
purities. Then, the inhomogeneity is expected only in the 
larger scale. 

We understand these results by considering the two 
typical scales of impurity concentration. One is the crit- 
ical concentration defined by the Abrikosov-Gorkov the- 
Qj.y 49-52 rpj^g other is defined so that the mean length 
between impurities becomes comparable with the coher- 
ence length in the clean limit, namely ^o- The former is 
scaled as nimp oc T^/Ef where Ef is the Fermi energy. 
The latter is obtained as riimp ~ oc (T^/Ef)'^ in two 
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Fig. 3. Spatial dependence of d-wave SC order parameter at T : 
. (a) The case of long coherence length superconductor, V - 



—0.6 and n^^p 
superconductor, 
^^imn = 0.05. 



= 0.01. (b) The case of short coherence length 
V = -1.5 and ni^p = 0.01. (c) V = -1.5 and 



dimension. The latter is much smaller than the former 
in the low temperature superconductors, and therefore 
the system has many impurities in the area ~ at the 
moderate impurity concentration. This is the underlying 
reason why the Abrikosov-Gorkov theory is justified. 

In contrast to the conventional case, the two scales 
can be comparable in the short coherence length super- 
conductor, such as the high-Tc and organic supercon- 
ductors. Then, the local superconductivity is stabilized 
in the "clean region" where only a few impurities ex- 
ist. For example, we show the distribution of point dis- 
orders for niinp = 0.05 in Fig. 4(a). Comparing with 
Fig. 3(c), we see that the SC order parameter is en- 
hanced in the nanoscale "clean region" which is shown 
by the red circles. The SC order parameter in Fig. 3(c) 
is localized in one of the red circles in Fig. 4(a). We un- 
derstand that the SC correlation develops in the other 
"clean regions" by showing A2(?) = i^) 
and A^li) = X](5(^l)^^(3i h ^) which are the d-wave or- 
der parameters having the second and third highest Tc, 
respectively. It is shown in Figs. 4(b) and (c) that they 
are developed in one of the red circles in Fig. 4(a). Thus, 
the highly disordered d-wave superconductor with short 
coherence length has a granular-like spatial structure. It 
should be stressed that the length scale of each gran- 
ular is nanometer scale which is much smaller than the 
conventional granular systems with macroscopic or meso- 
scopic scale. 

The microscopic inhomogeneity in Figs. 3(b), 3(c) and 
4 is consistent with the STM measurement in the under- 
doped region. This experimental observation is the fin- 
gerprint of the d-wave short coherence length supercon- 
ductor with point-like disorders. This interpretation is 
consistent with the recent experimental report^^ where 
the spatial structure of SC order parameter is associated 
with the extended but point impurities. Then, the clean 
region is characterized by the small gap and the sharp 
gap edge.^^ The broad gap structure in the dirty region 
should be attributed to the competing order. ^^'^^ 

As a result of the nanoscale inhomogeneity of SC order 
parameter, the mean field value of transition temperature 
is robust for the disorder, as pointed out in Ref. 58. 
This is simply because the localized order parameter in 
the "clean region" is little affected by the disorder poten- 
tial. We numerically take the random average of transi- 
tion temperature and show the result in Fig. 5. The 
result of SCTMA is also shown for a comparison. The de- 
crease of the mean field transition temperature remark- 
ably deviates from the SCTMA in case of the large \V\, 
namely the short coherence length {V = —1.5) while the 
deviation is negligible in case of the long coherence length 
{V — —0.6). Thus, the Abrikosov-Gorkov theory breaks 
down in the presence of the nanoscale inhomogeneity. 

Before closing this subsection, we note that the spa- 
tial dependence of SC order parameter is quite different 
in the s-wave superconductor. Then, the order parameter 
is almost uniform for 10% disorder-doping even if the co- 
herence length is short. This is because the s-wave su- 
perconductivity is robust for the non-magnetic impurity 
as proved by Anderson. ^° We obtain the spatial depen- 
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Fig. 4. (a) The distribution of point disorders in Fig. 3(c). We 
show the "clean region" by red circles. The SC correlation is lo- 
cally enhanced in these regions, (b) Spatial dependence of A2(j) 
and (c) that of A3(i). We choose the same parameters as in 
Fig. 3(c). 



dence similar to Fig. 3(a) where the SC order parameter 
is nearly homogeneous except for the abrupt suppression 
around impurities. Thus, the d-wave pairing symmetry 
plays an essential role in the inhomogeneous structure of 
SC order parameter in the nanometer scale. This may 
be the reason why the critical behavior in short coher- 
ence length s-wave superconductor NbN is quite different 
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Fig. 5. Random average of the transition temperature in 
the mean field theory. Here, we take more than 50 samples for 
the random average. The attractive interaction is varied as V = 
—0.6, —0.8, —1, —1.25 and —1.5 from the bottom to the top. 
We show the results of SCTMA with V = -0.6 and V = -1.5 
(dashed lines) for a comparison. 



from that in the under-doped cuprates.®^ 

3.2 SC fluctuation in the inhomogeneous system 

The purpose of this paper is to investigate the SC fluc- 
tuation in the disordered system with focus on the role 
of microscopic inhomogeneity. We discuss this issue by 
showing the results of self-consistent T-matrix approxi- 
mation. 

First, we discuss the transition temperature of super- 
conductivity. Because the long range order is not realized 
at finite temperature in purely two-dimensional systems 
according to the Mermin- Wagner theorem, we adopt a 
phenomenological procedure to include the weak three 
dimensionality. For this aim, the criterion |V^|ii = 1 — 5 
has been adopted instead of \V\ti = 1 when the clean 
system has been investigated.^'' We also adopt this pro- 
cedure and determine the Tc for nimp = by choosing 
6 = 0.02. However, it is not clear whether this criterion is 
applicable to the disordered systems. In particular, this 
procedure may break down in the disordered short coher- 
ence length superconductor with nanoscale inhomogene- 
ity. Then, it is expected that the long range correlation 
hardly develops although the local superconductivity oc- 
curs as in Figs. 3 and 4. Therefore, we determine the Tc 
by the following procedure. First, we calculate the cor- 
relation function of d-wave SC order parameter Xdi'P), 
which is defined as. 



Xd(*,j) = > (-ir+^mj,<5,<5'). 



Xd{f) =< 



N ^ 



Xd(j + >r, 



(22) 



(23) 



where Or denotes the random average. Here, P (P') 
is determined by j (j') as in eq. (20). We denote the 



correlation function for riimp = 
X;^ (f) and adopt the criterion, 

(c). 



and \Y\tx = 0.98 as 



Xd(15,15) = x:i'(15,15), 



(24) 
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to determine the in disordered systems. Note that 
f = (15, 15) is the longest length in the calculation for 
= 31 X 31 with (anti-)periodical boundary condition. 
This criterion determines the temperature where the long 
range correlation begins to develop. 




0.05 0.1 0.15 

n. 

imp 




0.05 0.1 0.15 

n. 

imp 

Fig. 6. The transition temperature obtained in the self-consistent 
T-matrix approximation for V = —1.5. (a) The results of the 
formulation in §2.1 where the disorder is exactly taken into ac- 
count. We take at least 20 samples and take the random average 
to estimate the correlation function, eq. (23). (b) The results of 
SCTMA. We show the transition temperature in the mean field 
theory, f^^ and T^^ by the dashed lines. 



We show the results of Tc in Fig. 6(a). The results ob- 
tained in the SCTMA are shown in Fig. 6(b) for a com- 
parison. The same criterion, namely eq. (24) is adopted 
in the SCTMA. The transition temperatures in the mean 
field theory ((a) f^'^ and (b) T^^ , respectively) are 
shown to illuminate the role of SC fluctuation. 

Since the short range SC fluctuation develops below 
the mean field transition temperature, can be re- 

garded as T* which is the onset of pseudogap phenom- 
ena arising from the SC fluctuation. An important value 
is the difference — which measures the width of 
critical region. When the SC fluctuation is strong, the 
critical region is large and vice versa. 

We clearly see the qualitatively different behaviors of 
critical region between Figs. 6(a) and (b). The results of 



the formulation in §2.1 show that the critical region is en- 
largedhy the disorder. However, the SCTMA provides an 
opposite result. This discrepancy means that the result of 
Fig. 6(b) is an artifact of the SCTMA because the disor- 
der is exactly taken into account in Fig. 6(a). This com- 
parison illuminates the breakdown of Abrikosov-Gorkov 
theory for the fluctuation phenomena in the short coher- 
ence superconductor.^^ 

The result of SCTMA, where the SC fluctuation is sup- 
pressed by the disorder, is seemingly reasonable because 
the Tc is decreased by the disorder. However, the micro- 
scopic inhomogeneity neglected in the SCTMA enhances 
the fluctuation and disturbs the long range order. The 
latter effect is significant in the short coherence length 
superconductor like under-doped cuprates and organic 
superconductor k-(BEDT-TTF)2X. Thus, the enhance- 
ment of SC fluctuation indicated by the experimental 
results"^""^^ is explained by appropriately taking into ac- 
count the inhomogcneous structure of SC fluctuation. 




Fig. 7. The spatially averaged DOS p{uj) for V = -1.5. (a) The 
results of the formulation in §2.1. We take more than 20 samples 
to take the random average, (b) The results of SCTMA. 



The experiments'^^"^^ have shown that the pseudogap 
is robust for the disorder even though the transition tem- 
perature is significantly decreased. In order to illuminate 
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this feature, we show the single particle DOS at T = T^. 
We calculate the local DOS p{r, uj) in each sample which 
is defined as, 

p{f, Lj) = --IniG^(f, f, w), (25) 

where the retarded Green function is obtained by the 
Fade approximation. By taking the spatial average and 
random average, we obtain the DOS per Cu sites as, 

pH=<-^^p(f,^)>r. (26) 
iVo — 

r 

As shown in Fig. 7(a), the pseudogap in the DOS is not 
significantly destroyed even though the is remark- 
ably reduced. This result should be contrasted to that 
of SCTMA where the pseudogap is smeared by the same 
amount of disorder (Fig. 7(b)). Thus, the robust pseu- 
dogap indicated by the experiments is attributed to the 
role of microscopic inhomogeneity. 




A| , 1 , 1 , 1 , 1 , 1 

1) 2 4 6 8 10 
X 

Fig. 8. The correlation function of SC order parameter Xdif)- 
We show the spatial dependence along the diagonal line f = 
{x,x). The solid line is the result of the calculation formulated 
in §2.1 and the dash-dotted line shows the result of SCTMA for 
"imp = 0.104. The dashed line is the result in the clean limit. 
We choose V = —1.5 and the temperature which satisfies the 
criterion eq. (24). Note that the correlation functions at large x 
are similar to each other owing to the criterion. 



In order to understand the results in Figs. 6 and 7 
more clearly, we again stress that the short range SC 
correlation develops even though the long range coher- 
ence is disturbed by the localization of order parame- 
ter. Because we have defined the as the temperature 
where the long range correlation begins to develop, the 
short range correlation at that temperature is enhanced 
by the disorder. This feature is clearly shown in Fig. 8 
which shows the correlation function of SC order parame- 
ter, namely Xdi'P)- The pseudogap in the short coherence 
length superconductor is robust for the disorder since the 
pseudogap is mostly induced by the short range correla- 
tion.^'' We note that the qualitatively different results 
are obtained in the SCTMA where the short range cor- 
relation is reduced by the disorder. This is because the 
coherence length increases owing to the decrease of Tc 
and the long range correlation tends to develop. This 
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feature is consistent with Fig. 7(b), but an artifact of 
the SCTMA. 

It should be noted that the long range correlation is 
still weak in Fig. 8. This means that the actual long range 
coherence is achieved at the temperature lower than Tc 
in our definition. In order to discuss the true long range 
order, we have to take the weak three dimensionality into 
account. The nature of dimensional crossover from 2D to 
3D has been investigated in the clean systems. Then, the 
qualitative behaviors of Tc do not depend on the phe- 
nomenological criterion.*^ However, the nature of long 
range order is not clear in the highly disordered system. 
It is expected from Figs. 3 and 4 that the long range order 
is triggered by the phase coherence between each "clean 
region" like granular superconductivity.^^ The result in 
Fig. 8 is consistent with this picture although the calcu- 
lation with = 31 X 31 is not sufficient to describe the 
phase correlation between each granular. We have con- 
firmed that the behaviors in Figs. 6 and 8 become notable 
as increasing the size of calculation. It is an interesting 
future problem to develop a theory on the crossover from 
the homogeneous superconductivity to the granular su- 
perconductivity. 

4. Pseudogap around Single Impurity 

In §3, we have investigated the role of microscopic in- 
homogeneity on the macroscopic properties of d-wave 
short coherence length superconductor. We briefiy dis- 
cuss the spatial dependence in the atomic scale in this 
section. 

Here, we focus on the spatial dependence of the elec- 
tronic DOS in the presence of strong point disorders like 
Zn impurities. Because the most significant spatial de- 
pendence is the variation around the disorder, we show 
the DOS around the single impurity. The pseudogap 
around the single impurity has been investigated by the 
NMR measurements. ^''"^^ We discuss the relevance of our 
results for these experiments. 



0.5 




Fig. 9. The LDOS around the single impurity. We show the 
LDOS at the nearest neighbor (NN) site (solid line) and next 
nearest neighbor (NNN) site (dot-dashed line) in addition to the 
spatially averaged DOS (dashed line). 

In this section, we adopt the self-consistent T-matrix 
approximation for = 31 x 31 sites. Figure 9 shows the 
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local DOS (LDOS) at the nearest neighbor (NN) and 
next nearest neighbor (NNN) sites where the deviation 
from the averaged DOS is most significant. The temper- 
ature is chosen to be slightly above Tc and therefore the 
pseudogap opens in the averaged DOS as shown by the 
dashed line. It is clearly shown that the pseudogap is de- 
stroyed at the NN site. This is because the Cooper pair- 
ing is forbidden between the impurity site and NN site. 
We have confirmed that the d-wave order parameter is re- 
markably suppressed at the NN site. The destroyed pseu- 
dogap at the NN site is consistent with the NMR mea- 
surements which have shown the Curie- Weiss enhance- 
ment of l/TiT around the non-magnetic impurity.^'*"^^ 
The neutron scattering measurements have also shown 
the destruction of superconducting gap in the magnetic 
excitation.^*' 
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Fig. 10. The spatial dependence of the LDOS along the [100] and 
[110] directions from the impurity site. We define the impurity 
site f = (0,0) and show the LDOS p{f,0) at r = {x,0) (circles) 
and f = {x,x) (squares), respectively. 



Figure 9 shows that the pseudogap in the DOS at the 
NNN site is larger than that in the averaged DOS. Thus, 
the spatial dependence of the DOS is non-monotonic. We 
show the LDOS p(f, 0) at the Fermi energy along [100] 
and [110] directions in Fig. 10. The DOS remarkably os- 
cillates along [100] direction while the DOS shows a long 
wave length oscillation along [110] direction. It should 
be noted that the SC order parameter is monotonically 
suppressed around the impurity site and shows no os- 
cillation. Thus, the suppression of the DOS along the 
[110] direction is not attributed to the enhancement of 
the SC correlation. Actually, the oscillation of the DOS 
is a kind of the quasi-particle interference effect^° which 
reflects the shape of the Fermi surface. 

We have analyzed the magnetic excitations in the dis- 
ordered t-t'-U-V model in order to take into account 
the competing anti-ferromagnetic fluctuation^^ in ad- 
dition to the SC fluctuation. Then, we found that the 
spatial dependences of static quantities are quite differ- 
ent from that of NMR 1/TiT. The staggard spin sus- 
ceptibility is monotonically enhanced near the impurity. 
This enhancement of the anti-ferromagnetic correlation 
induces the rapid oscillation of uniform spin susceptibil- 



ity which has been observed in the Knight shift mea- 
surements.^^ This spatial dependence of the static anti- 
ferromagnetic correlation can be qualitatively described 
by the Ginzburg-Landau theory which assumes the com- 
petition between the anti-ferromagnetism and supercon- 
ductivity. In contrast to that, the \/TiT shows a spa- 
tial dependence similar to the electronic DOS in Fig. 10, 
which can not be described in the Ginzburg-Landau the- 
ory. This is consistent with the experimental result on 
the l/TiT which shows the length scale shorter than the 
Ginzburg-Landau correlation length. ^"^ The details will 
be shown in the future publication.''^^ 

When we distribute many strong point disorders like 
Figs. 3 and 4, the Fourier transformed DOS (FTDOS) is 
dominated by the deviation at the NN and NNN sites. 
Then, no clear structure appears in the FTDOS around 
q — (7r/4, 0). In contrast to that, the peak appears in the 
FTDOS around q — (7r/4, 0) (not shown) for the weak 
random disorder. The latter is similar to the experimen- 
tal observation in the pseudogap state, although the 
structure around q — (7r/4, 0) is smeared in the high fre- 
quency region. 

We note that the nanoscale inhomogeneity of DOS 
observed in the STJVI^^"^^ is not clearly reproduced in 
this calculation. We consider that the huge and broad 
gap structure in the weakly superconducting region is 
due to the correlation competing with the d-wave su- 
perconductivity.^^ A candidate is the disordered mag- 
netism observed in the neutron, /iSR and N]VIR mea- 
surements. ^^"^^ 

5. Discussion 

5.1 Summary and discussion 

We have investigated the disorder-induced microscopic 
inhomogeneity of SC order parameter and its effects on 
the pseudogap in short coherence length superconductors 
with d-wave symmetry. We find that the SC order param- 
eter is highly localized by a small amount of point dis- 
orders and the nanoscale inhomogeneity like a granular 
superconductor appears. Owing to the spatially localized 
SC correlation, the long range coherence hardly develops 
while the short range correlation is enhanced. As a re- 
sult, the critical region of superconductivity is enlarged 
by the disorder in contrast to the SCTJVIA based on the 
Abrikosov-Gorkov theory.^*' Then, the pseudogap in- 
duced by the SC fluctuation is robust for the disorders. 

One of the typical short coherence length supercon- 
ductors is the high-Tc superconductor. The under-doped 
cuprates satisfy three conditions to realize the micro- 
scopic inhomogeneity, namely the short coherence length, 
non-s-wave symmetry and quasi-two dimensionality. Ac- 
tually, the nanoscale inhomogeneity has been observed in 
the STIVE measurements^^"^^ and also in the Josephson 
plasma resonance measurement.^^ We note that the dis- 
order induces the inhomogeneity of SC order parameter 
without remarkable charge inhomogeneity. The results 
of NIVIR and transport measurements^^"^^'^^ indicating 
the breakdown of Abrikosov-Gorkov theory are consis- 
tent with the microscopic inhomogeneity induced by the 
disorder. We understand these experimental results in a 
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coherent way by appropriately taking the spatial depen- 
dence into account. 

It should be stressed that there are intrinsic disorders 
in high-Tc cuprates. Therefore, our results indicate the 
nanoscale inhomogeneity in the non-disorder-doped sam- 
ple. Then, it is expected that the superconducting tran- 
sition in the highly under-doped region is triggered by 
the phase coherence of nanoscale granular regions. This 
picture is consistent with the Uemura plot''^ which has 
implied the importance of phase fluctuation as pointed 
out by Emery and Kievelson.^ The recent experiments 
have also indicated the two dimensional phase fluctua- 
tion above T^.^'^''^'^''^^ 

The nature of the SI transition in the highly under- 
doped region is beyond the scope of this paper. How- 
ever, it is reasonable to consider that the disordered mag- 
netism occurs as shown by the experiments. ^^"^^ We will 
discuss this topic in another publication.^^ 

5.2 Discussion on a variety of SI transitions 

A variety of SI transitions have been observed in 
strongly correlated electron systems. Here, we briefly dis- 
cuss how the differences can be midcrstood by consider- 
ing the role of disorder. The importance of disorder in 
the under-doped cuprates is illuminated by the following 
discussion. 

The nature of SI transition is different between the 
series of high-Tc cuprates. The disordered magnc;tism 
has been observed most clearly in the La-based systems 
where the Cu02 layer is significantly affected by the dis- 
order. In contrast to that, the recent work on the multi- 
layer cuprates has shown the abrupt change from the 
superconducting state to the anti-ferromagnetic metal- 
lic state at finite temperatures.^® The latter is consis- 
tent with the theoretical results based on the FLEX and 
FLEX-|-T-matrix approximations in the clean system. 
The authors have claimed that the clean Cu02 layer is 
realized in this compound. ^® Then, the absence of the 
disordered magnetism and/or the anti- ferromagnetic in- 
sulating state implies that the disorder plays an essential 
role for the SI transitions in the other cuprates. 

The organic superconductor k-(BEDT-TTF)2X*° is 
another c;andidatc for the short coherence length d-wave 
superconductor.^^' It is expected that the system is rel- 
atively clean because the material is free from the dopant 
disorders. The first order phase transition from the 
superconductivity to the anti-ferromagnetism indicates 
that the organic superconductor k-(BEDT-TTF)2X is 
analogous with the multi-layer high-Tc cuprates. 

The electron-hole asymmetry of the SI transition in 
cuprate superconductors is understood by considering 
the difference of coherence length. The result of the 
FLEX-l-T-matrix approximation has shown that the co- 
herence length is much longer in the electron-doped re- 
gion than in the hole-doped region.'''* This is the rea- 
son why the pseudogap due to the SC fluctuation is not 
observed in the electron-doped cuprates. ^^"^^ Then, the 
abrupt change from the superconducting state to the 
anti-ferromagnetic state*^ is attributed to the absence of 
microscopic inhomogeneity which is unlikely in the long 
coherence length superconductor. 



A variety of SI transitions accompanying the mag- 
netism have been discovered in the heavy fermion mate- 
rials.^^ It is reasonable to consider that the multi-orbital 
effect plays an essential role for the phase transition from 
the superconductivity to the magnetism. In some cases, 
the superconductivity and magnetism are induced by the 
different orbitals. Then, the coupling between the super- 
conductivity and magnetism is weak and the co-existent 
state can be stabilized. 

5. 3 Plan of the future study 

It is a challenging future work to describe the SI transi- 
tion on the basis of the disorder-induced nanoscale inho- 
mogeneity.*^ Some theoretical developments are needed 
for this aim. 

First, it is important to take into account the 
competition between the superconductivity and anti- 
ferromagnetism. This competition induces the spatially 
inhomogeneous structure of superconducting region and 
magnetic region.®* We have examined the spatial struc- 
ture on the basis of the t-t'-U-V model which includes 
the on-site repulsion. Then, we found that the magnetic 
correlation is enhanced in the "dirty region" where the 
superconducting correlation is destroyed by the disorder. 
The results will be shown in the future publication.^^ 

Second, it is necessary to take into account the quan- 
tum fluctuation to describe the SI transition at zero tem- 
perature although we have investigated the thermal fluc- 
tuation in this paper. One of the promising strategies 
is to derive the phenomenological models such as the 
Ginzburg-Landau model or the XY model where the long 
range behaviors are illuminated.*® We are planning to 
microscopically derive a phenomenological model which 
is relevant in the spatially coexistent state of the super- 
conductivity and anti-ferromagnetism. 
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